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Abstract 

The pion wave function is discussed in the hght of the recent CLEO 
data on the 7r7 transition form factor. It turns out that the wave function is 
close to the asymptotic form whereas wave functions strongly concentrated 
in the end-point regions are disfavoured. Consequences for other exclusive 
quantities, as for instance the pion's electromagnetic form factor, are also 
discussed. 
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The theoretical description of large momentum transfer exclusive reactions is 
based on a factorization of long- and short- distance physics. The latter physics 
is contained in the so-called hard scattering amplitude to be calculated within 
perturbation theory. Universal, process-independent (light cone) wave functions 
interpolating between hadronic and partonic degrees of freedom, comprise the 
long-distance physics. The wave functions are not calculable with sufficient de- 
gree of accuracy at present. However, for the pion which is the hadron of interest 
in this letter, the valence Fock state wave function is theoretically rather well con- 
strained. The main uncertainty of the wave function lies in the x-dependent part 
of it, the so-called distribution amplitude x denotes the usual momentum 

fraction the valence quark carries. 

Previous studies of the pion's electromagnetic form factor as well as other 
large momentum transfer exclusive reactions, as for instance 77 — vrvr, seemed to 
indicate that the pion distribution amplitude is much broader than the so-called 
asymptotic one, ipAsi^) = 6x(l — x). Chernyak and Zhitnitsky [|I| proposed 
such a broad distribution amplitude which is strongly end-point concentrated 
and leads to a leading twist contribution to the pion's electromagnetic form factor 
in apparently fair agreement with the admittedly poor data 0. This result is, 
however, obtained at the expense of the dominance of contributions from the end- 
point regions, x or 1, where the use of perturbative QCD is unjustified as has 
been pointed out by several authors ^. Now the prevailing opinion is that the 
pion's electromagnetic form factor is controlled by soft physics (e.g. the overlap 
of the initial and final state wave functions, occasionally termed the Feynman 
contribution) for momentum transfer less than about 10 GeV^ ||^, |^, |^, 0]. 

There is another exclusive quantity namely the 717 transition form factor 
which, for experimental and theoretical reasons, allows a more severe test of 
our knowledge of the pion wave function than the pion's electromagnetic form 
factor. Recently the form factor has been measured in the momentum transfer 
region from 2 to 8 GeV^ with rather high precision. Together with previous 
CELLO measurement [0| we now have at our disposal much better data above 1 
GeV^ for the irj transition form factor than for the pion form factor. From the 
theoretical point of view the analysis of the 717 transition form factor is much 
simpler than that of the pion form factor: It is, to lowest order, a QED process, 
QCD only provides corrections of the order of 10% in the momentum transfer 
region of interest |10[ . The difficulties with the end-point regions where the gluon 



virtuality becomes small, do not occur. Higher Fock state contributions, sup- 
pressed by powers of a^/Q^, are expected to be small (see the discussion in |ri[]). 
Moreover, and in contrast to the pion form factor, the Feynman contribution, 
which may arise through vector meson dominance, is presumably very small due 
to a helicity mismatch. 

The purpose of this letter is to extract information on the pion wave function 
from a perturbative analysis of the 717 form factor. It will turn out that the new 
CLEO data p allow a fairly precise determination of that wave function. It will 
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also be argued that this wave function leads to a consistent description of the 
pion's electromagnetic form factor and its structure function. It should be noted 
that this letter is an update of previous work [|T^. The very important CLEO 
data H] were not yet available in [O. We also point out that the tt'j* transition 



form factor is investigated in |]13 



Let us begin with the parameterization of the soft valence Fock state wave 
function, i.e. the full wave function with the perturbative tail removed from it. 
This is the object required in a perturbative calculation [0. Following the 
soft wave function is written as 



^0 {x,h,fiF) 



2V6 



X, jj^p)^ ( y x{l — x)h 



(1) 



is the scale at which soft and hard physics factorize |T^, |T5|, [T^. The wave 
function is subject to the auxiliary conditions 



S(0) = in, 



dx(f){x, fip) = 1. 



(2) 



b is the quark-antiquark separation in the transverse configuration space and is 
canonically conjugated to the usual transverse momentum k^. It is advanta- 
geous to work in the transverse configuration space because the Sudakov factor, 
to be discussed later, is only derived in that space (see [l^). The parameter- 
ization (0) automatically satisfies the constraint from the process n^ f^^^^i 
Q which relates the wave function at the 'origin of the configuration space' to 
the pion decay constant /7r(= 130.7 MeV). On the assumption of duality prop- 
erties, one can derive constraints on x and k± moments^ of the wave function 
(e.g. (x") = J dxx"^ 4'{x)) within the operator product expansion framework ]T7|. 
These constraints can be combined into the following conditions on the momen- 
tum space wave function: 

i) The distribution amplitude has simple zeroes at a; 0, 1. 

ii) The k^-dependence of the wave function \l/o {x, k^) comes exclusively 
in the combination k'j_/x{l — x) at a; 0, 1. 

iii) At large k± the wave function falls off faster than any power of k^. 



The simplest function matching the conditions ii) and iii) is the Gaussian 

x(l - x) 6^ 



x(l — x)h 



47rexp 



4a2 



(3) 



where a is the transverse size parameter. This Gaussian will be used subsequently. 



^ There is a subtlety: The kinematical transverse momentum of the parton is not the same 
object as k± defined through the moments. In this letter we will assume that both are one and 
the same variable. This assumption corresponds to summing up soft gluon corrections, i.e. to 
higher twist contributions. 
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The distribution amplitude is subject to evolution and can be expanded over 
Gegenbauer polynomials C^''^, the eigenfunctions of the (leading order) evolution 
equation for mesons 



1. t 



?i=2,4,... 



(4) 



«s is the strong coupling constant and /iq is a typical hadronic scale for which we 
choose 0.5 GeV throughout. Charge conjugation invariance requires the odd n 
expansion coefficients Bn to vanish. Since the 7„ are positive fractional numbers 
increasing with n (e.g. 72 = 50/81) any distribution amplitude evolves into (pAs 
asymptotically, i.e. for ln((5//io) 00; higher order terms in (|) are gradually 
suppressed. The asymptotic distribution amplitude itself shows no evolution. 
This property of (pAs no more holds if evolution is treated in next-to-leading or- 
der As we are going to show below the asymptotic distribution amplitude 



in combination with the Gaussian (|^) provides very good results for the 717 form 
factor and also leads to a consistent description of other exclusive reactions in- 
volving pions. In order to quantify the amount of deviations from the asymptotic 
distribution amplitude still allowed by the CLEO data, contributions from the 
second Gegenbauer polynomial will be permitted in the analysis and limits for the 
strength of the coefficient B2 will be extracted. For the purpose of comparison the 
7f~f form factor will be evaluated with the Chernyak-Zhitnitsky (CZ) distribution 
amplitude IQ which is defined by i?2 = 2/3 and 5„ = for n > 2. 

For a given distribution amplitude there is only one free parameter in the 
wave function, namely the transverse size parameter a. It can be fixed by using 
a constraint derived from the process —>■ 77 [0]: 



dx(fh'^o{x,h) = —. (5) 

Jtt 

Although not at the same level of rigour as the other constraints (because of 
approximations made in the case when one photon couples 'inside' the pion wave 
function) (^) still provides a value for the parameter a which comes up to our 
expectations for the transverse size of the pion. In particular for the asymptotic 
wave function one finds a value of 861 MeV for a, which corresponds to a value 
of 367 MeV for the root mean square transverse momentum. The probability of 
the valence Fock state amounts to 0.25 in this case. 

Now the wave function is fully specified and we can turn to the calculation of 
the 7r7 form factor. Inspection of the ttj form factor data [^, Q (see Fig. 1) reveals 
a dependence somewhat stronger than predicted by dimensional counting. 
Hence, higher twist contributions do not seem to be negligible in the momentum 
transfer region from 1 to 8 GeV^. The modified perturbative approach proposed 



by Sterman and collaborators [|lg, |T9| allows to calculate some power corrections 
to the leading twist term. In that approach the transverse momentum depen- 
dence of the hard scattering amplitude is retained and Sudakov suppressions are 
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taken into account in contrast to the standard approach Apphcations of 

the modified perturbative approach to the pion's and nucleon's electromagnetic 
form factors ^ revealed that the perturbative contributions to these form 

factors can self-consistently (in the sense that the bulk of the contributions is 
accumulated in regions where the strong coupling is sufficiently small) be cal- 
culated. It turned out, however, that the perturbative contributions are too small 
as compared with the data. Therefore, in the experimentally accessible range of 
momentum transfer, these form factors are controlled by soft physics. Higher 
order perturbative corrections and/or higher Fock state contributions seem too 
small in order to account for the large discrepancies between the lowest order 



perturbative contributions and the data for the elastic form factors (see |^ for 
the pion case). 

For the reasons discussed in the introduction the 717 transition form factor 
represents an exceptional case for which we can expect perturbation theory to 
work for larger than about 1 or 2 GeV^. Adapting the modified perturbative 
approach |]16|, |l^ to the case of transitions, we write the corresponding form 



factor as 

iff) = J dx {x, -b, jj^p) Th (x, b, Q) exp [-S (x, 6, Q)] . (6) 



This convolution formula can formally be derived using the methods described 
in detail by Botts and Sterman DTG). Th is the Fourier transform of the momen- 
tum space hard scattering amplitude to be calculated, to lowest order, from the 
Feynman graphs shown in Fig. 2. It reads 

Th (x, b, Q) = -^Ko (VT^ Q b) (7) 

where Kq is the modified Bessel function of order zero. The Sudakov exponent 
in (^, comprising those gluonic radiative corrections (in next-to-leading-log 
approximation) not taken into account in the evolution of the wave function, is 
given by 

S{x, b, Q) = s{x, b, Q) + s(l - X, 6, g) - ^ In ■ (^/^^^^^ 



(3o ln(l/6AQCD) 



where a Sudakov function s appears for each quark line entering the hard scat- 
tering amplitude. The last term in (H) arises from the application of the renor- 
malization group equation {[3q = 11 — |n/). A value of 200 MeV for Aqcd is used 
and fi is taken to be the largest mass scale appearing in the hard scattering am- 



plitude, i.e. /i = max — xQ, 1/bj. For small b there is no suppression from 
the Sudakov factor; as b increases the Sudakov factor decreases, reaching zero at 
b = 1/Aqcd- For even larger b the Sudakov is set to zero. The Sudakov function 



s has been calculated by Botts and Sterman [|T6] using resummation techniques 
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its explicit form can be found in Due to the properties of the Sudakov fac- 
tor any contribution is damped asymptotically, i.e. for ln((5^//io) oo, except 
those from configurations with small quark-antiquark separations and, as can be 
shown, the hmiting behaviour F^^^ — > V^fw/Q"^ emerges, a result which as been 
derived previously |]23| , |2^ . b plays the role of an infrared cut-off; it sets up the 
interface between non-perturbative soft gluon contributions - still contained in 
the hadronic wave function - and perturbative soft gluon contributions accounted 
for by the Sudakov factor. Hence, the factorization scale /li? is to be taken as 1/6. 

Employing the asymptotic wave function, (|l|), (|) and (j)As, "we find, from 
(^, the numerical results for the 7r7 transition form factor displayed in Fig. 1. 
Obviously there is very good agreement with the data || above — 1 GeV^. 
At 8 GeV^ about 85% of the asymptotic value has been reached. We emphasize 
that there is no free parameter in our approach to be fitted to the data once the 
wave function is chosen and the transverse size parameter is fixed through (^. For 
comparison we also show in Fig. 1 results obtained with the CZ wave function, (P, 
(H) and (pcz (-B2 = 2/3, = for n > 2 in (H)). That prediction overshoots the 
data markedly. Of course, the experimental errors allow slight modifications of the 
asymptotic wave function. In order to give a quantitative estimate of the allowed 
modifications we fit the expansion coefficient B2 to the data assuming 5„ = for 
n > 2 and choosing again fiQ = 0.5 GeV. For each value of B2 the transverse size 
parameter a is fixed through (^. A best fit to the data above 1 GeV^ provides 
^mpa ^ _o.o06 ± 0.014 (with a = 864 MeV), i.e. a value compatible with 
zero. In [|14| a modification of the asymptotic wave function is proposed where 
(pAs is multiplied by the exponential exp — m^a^/x(l — x) . The parameter rriq 
represents a constituent quark mass of, say, 330 MeV. A similar wave function 
is constructed by Dorokhov from the helicity and fiavour changing instanton 
force. Although wave functions of this type contradict the constraint i) derived 
by Chibisov and Zhitnitsky [|17|, they cannot be excluded absolutely since the 



constraint i) is obtained under a duality assumption the validity of which is not 
guaranteed. In any case we have convinced ourselves that the wave function given 
H[ provides similarly good results for the 717 form factor as the asymptotic 



m 



wave function itself. This is not a surprise since both the wave functions differ 
from each other only in the end-point regions. Contributions from these regions 
are strongly suppressed by the Sudakov factor. 

It is instructive to compare the leading twist result for the 717 form factor 



F., {Q') = ^ {X-') 1^ [ 1 + ^^K{Q, +0{al)] (9) 

with the data above ^ 3 GeV^ which, within errors, are just compatible with 
a Q^-dependence according to dimensional counting, /i/j represents the renormal- 
ization scale. The factor K depends on the distribution amplitude. Using the 
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expansion (^), one finds for tlie x ^ moment of the distribution amplitude 



(10) 



Neglecting again terms with n > 2 and also the corrections in (^), we obtain 
a value of —0.39 ± 0.05 for -B^''^ from a fit to the data {fip = Q in this case). The 
face value of -B^''^ corresponds to {x~^)^^ = 2.39 (at Q'^ = 8 GeV^) which is to 
be contrasted with the values of 3 and 4.01 (at 8 GeV^) for the asymptotic and 
the CZ distribution amplitude respectively. 



Braaten |T^ has calculated the corrections (in the MS scheme) in (^). His 
analysis is however incomplete in so far as only the corrections to the hard 
scattering amplitude have been considered but the corresponding corrections to 
the kernel of the evolution equation for the pion's distribution amplitude were 
ignored. As has been shown by Miiller |TB[ recently in next-to-leading order 



the evolution provides logarithmic modifications in the end-point regions for any 
distribution amplitude, i. e. for the asymptotic one too. An estimate however 
reveals that the modifications of the evolution behaviour in next-lo- leading order 
are very small for the asymptotic distribution amplitude (a^ evaluated in two-loop 
with A*^'?^3 = 200 MeV), and can safely be neglected here. For the CZ distribution 
amplitude these effects seem to be somewhat larger than for the asymptotic one 
but still the total as corrections are dominated by those to the hard scattering 
amplitude. The as corrections amount to —10% in the case of the CZ distribution 
amplitude. Hence, also in the leading twist analysis in next-to-leading order the 
CZ distribution amplitude is clearly at variance with the data. Next we want to 
determine the expansion coefficient B2 in the next-to-leading order leading twist 
analysis in order to quantify the deviations from the asymptotic distribution 
amplitude required by the data. For this purpose we evaluate the K factor 
in (H) from the expressions for the corrections given in [1^ (/i/j = jJ-F = Q) 
and neglect the modifications of the evolutions in next-to-leading order. From a 
fit of (^) to the form factor data we find B^^"-^ = —0.17 ± 0.05 corresponding 
to {x~^)^^^ = 2.74 at = 8 GeV^. According to what we said above such 
small a value of B2 will not be altered substantially by the modifications of the 
evolution behaviour to that order. Thus, the leading twist analysis requires a 
distribution amplitude which is a little narrower than the asymptotic one. In the 
modified perturbative approach, on the other hand, the asymptotic wave function 
works well since the QCD corrections condensed in the Sudakov factor, and the 
transverse degrees of freedom provide the required Q^-dependent suppressions. It 
is to be stressed that the Sudakov factor already takes into account the leading 
and next-to-leading logs of the as corrections. 

Other models, applicable at large as well as at low Q^, provide a parameteri- 
zation of the form factor as 



F^,(Q')=A/{l + Q'/so). (11) 
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Thus, Brodsky and Lepage |]23| propose that parameterization as an interpolation 
between the two hmits, Fj^-y (Q^ = 0) = A = (2A/27r^/,r)~^ known from current 
algebra and the limiting behaviour ^/2/^Q~^. Hence, sq = 47r^/^ = 0.67 GeV^ in 
that model. The interpolation formula works rather nicely, its Q^- dependence is 
similar to that one predicted by the modified perturbative approach. The vector 
meson dominance model leads to ([TT|) under the neglect of contributions from the 
meson and by ignoring the small mass difference between the p and u meson. 
The constant A is related to 

A=^ + ^ (12) 

fp fuj 

in the vector meson dominance model. Sq equals where rrip is the p meson 
mass. Inserting the known values of coupling constants [^, one finds for A a value 
of 0.269 ± 0.019 in agreement with the current algebra value. The vector meson 
dominance model is in accord with the present data although its asymptotic value 
{Arn^ ~ 0.16) differs from our one. A QCD sum rule analysis ||2^ also provides 



results similar to (|lTD . In order to discriminate among the various models data 
extending to larger values of momentum transfer and/or with smaller errors as 
the present ones are needed. 

Let us now turn to the discussion of the implications of our findings, namely 
that the perturbative analysis of the 7r7 transition form factor requires the asymp- 
totic pion wave function and apparently excludes strongly end-point concentrated 
wave functions like the one proposed by Chernyak and Zhitnitsky. Since the wave 
functions are universal, process-independent objects they should also be used in 
other large momentum transfer exclusive reactions involving pions, as for exam- 
ple the electromagnetic form factor of the pion or 77 tttt. As is well-known 
the leading twist results are only in agreement with experiment provided an end- 
point concentrated wave function respectively distribution amplitude is utilized^. 
This apparent agreement with experiment is, as we already mentioned, only ob- 
tained at the expense of strong contributions from the soft end-point regions 
where the use of perturbation theory is unjustified. This is to be contrasted with 
the modified perturbative approach where the end-point regions are strongly sup- 
pressed and a theoretically self-consistent perturbative contribution is obtained. 
However, as shown in ^ for the case of the pion form factor, the perturbative 
contributions evaluated with both the wave functions, the asymptotic one and 
the CZ one, are too small as compared with the data. At this point we remind 
the reader of the fact that the pion form factor also gets contributions from the 
overlap of the initial and final state soft wave functions \l/o ([1|). Formally the 
perturbative contribution to the pion form factor represents the overlap of the 
large momentum tails of the wave functions while the overlap of the soft parts of 



^ Taking the moment (x^^) which represents a soft, process-independent parameter, from 
our leading twist analysis of -FV7, we find a value for the pion form factor too small as compared 
to the admittedly poor data (note: F.^ ^ {x~^)^) even when corrections are considered pl|. 
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the wave functions is customarily assumed to be negligible at large Q. Examining 
the validity of that presumption by estimating the Feynman contribution from 
the asymptotic wave function, one finds results of appropriate magnitude to fill 
in the gap between the perturbative contribution and the data of Ref. 0]. The 
results exhibit a broad flat maximum which, for momentum transfers between 3 
and about 15 GeV^, simulates the dimensional counting behaviour]^. For the CZ 
wave function, on the other hand, the Feynman contribution exceeds the data 
significantly. Similar large Feynman contributions have also been obtained by 
other authors [@, 0]. Thus, the small size of the perturbative contribution to 
the elastic form factor finds a comforting although mo del- dependent explanation, 
a fact which has been pointed out by Isgur and Llewellyn Smith long time 
ago. 

The structure function of the pion offers another possibility to test the wave 
function against data. As has been shown in [|I^ the parton distribution func- 
tions are determined by the Fock state wave functions. Since each Fock state 
contributes through the modulus squared of its wave function integrated over 
transverse momenta up to Q and over all fractions x except those pertaining 
to the type of parton considered, the contribution from the valence Fock state 
should not exceed the data of the valence quark structure function. As discussed 
1^, ^ the asymptotic wave function respects this inequality while the CZ one 



m 



again fails dramatically. 

To conclude the asymptotic pion wave function, respecting all theoretical con- 
straints, provides a consistent and theoretically satisfying description of the nj 
and the pion's electromagnetic form factor and is compatible with the pion's 
valence quark distribution function. The pion's electromagnetic form factor is 
controlled by soft physics (which can be modelled as the Feynman contribution 
for the asymptotic wave function) in the experimentally accessible range of mo- 
mentum transfer in contrast to the transition form factor which is dominated 
by the perturbative contribution. We note that similar observations about the 
smallness of the perturbative contributions and the dominance of the Feynman 
contributions have been made in the case of the nucleon's form factor p8[. Of 
course, the present quality of the data does not allow to pin down the form of 
the wave function exactly. Little modifications of the asymptotic wave function 
can not be excluded as yet. On the other hand, the CZ wave function is in 
conflict with the data and ought to be discarded. This is also true for other 
strongly end-point concentrated wave functions. The use of such wave functions 
in the analyses of other exclusive reactions involving pions, e.g. 77 — > vrvr or 
B TTTT, seems to be unjustified (if one accepts the process-independence of the 
wave function) and likely leads to overestimates of the perturbative contributions. 



^ At large the Feynman contribution is suppressed by as compared to the pertur- 

bative contribution. The latter dominates the elastic form factor only for > 50 GeV^ . This 
value of is, however, very sensitive to the end-point behaviour of the wave function, little 
modifications may change it considerably. 
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The next-to-leading order leading twist analysis of the F^^^.^ form factor (possible 
for <?> GeV^) also reveals that the wave function or better the distribution 
amplitude in that case, is close to the asymptotic one but a little narrower than 
it. This result is to be contrasted with the modified perturbative approach where 
the asymptotic wave function works well; the required Q^-dependent suppression 
is provided by the Sudakov factor and the transverse degrees of freedom. In any 
case a systematic next-to-leading order analysis of exclusive reactions involving 
pions is required. 
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Figure captions 



Fig. 1. The scaled transition form factor vs. Q^. The sohd (dashed) hne 
represents the results obtained with the modified perturbative approach using the 
asymptotic (CZ) wave function. The evolution of the CZ wave function is taken 
into account. The dotted line represents the limiting behaviour \/2/vr. Data are 
taken from ||. 

Fig. 2. The basic graphs for the ttj transition form factor. 
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Fig. 2 



